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Abstract. The calculation of the conductance of disordered rings requires a 
theory that goes beyond the Kubo-Drude formulation. Assuming "mesoscopic" 
circumstances the analysis of the electro-driven transitions show similarities with a 
percolation problem in energy space. We argue that the texture and the sparsity 
of the perturbation matrix dictate the value of the conductance, and study its 
dependence on the disorder strength, ranging from the ballistic to the Anderson 
localization regime. An improved sparse random matrix model is introduced 
to captures the essential ingredients of the problem, and leads to a generalized 
variable range hopping picture. 



1. Introduction 

Closed mesoscopic rings provide the ideal paradigm for testing the manifestation of 
quantum mechanical effects in the mesoscopic realm [1] [2l [3l [4] . First measurements of 
the conductance of closed rings have been reported more than a decade ago [5j , while 
more recently there is a renewed experimental interest motivated by high precision 
measurements of individual rings 6, 7J. In a typical experiment a collection of 
mesoscopic rings is driven by a time dependent magnetic flux which creates an electro- 
motive-force (EMF). In what follows we assume low frequency DC noisy driving 
(u ~ 0) with power spectrum 



2u> c \ uj,_ _ 

where e is the RMS value of the voltage, and the cutoff frequency uj c is small compared 
with any relevant semiclassical energy scale, but larger compared with the mean level 
spacing A (note!). Optionally, if we had assumed an interaction with a thermal bath, 
the role of oj c would have been played by the level broadening or by the temperature 
[8]. In such setup one expects the rate of energy absorption to be given by Joule's law 
Ge 2 , where the coefficient G is defined as the "conductance'^. 

As in the linear response theory (LRT) analysis, we assume that the coherence 
time is much longer compared with the ballistic time, but smaller compared with the 

\ Hence there is no issue of quantum recurrences which would arise for a strictly linear or periodic 
driving. From here on we use units such that h = 1 . 

§ The terminology of this paper, and in particular our notion of "conductance" are the same as in 
the theoretical review [4] and in the experimental work [5|. 
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Heisenberg time (1/A). But our interest is in what we call mesoscopic circumstances. 
Namely, we assume that the environmental induced relaxation is a slow process, when 
compared with the EMF-driven transitions [this is in opposite to the LRT limit s — > 0]. 
Accordingly we have to work within the framework of semi-linear response theory 
(SLRT) [9l [lOl [TTfjj l: This theory (see section 5 below) goes beyond the conventional 
framework of the Kubo formalism. 

For diffusive rings the Kubo formalism leads to the Drude formula for G. A 
major challenge in past studies was to calculate the weak localization corrections [3] 
to the Drude result, taking into account the level statistics and the type of occupation 
[4] . These corrections are of order A/cu c and accordingly do not challenge the leading 
order Kubo-Drude result. It is just natural to ask what happens to the Drude result 
if the disorder becomes weak (ballistic case) or strong (Anderson localization case). 
In the latter case there are two conflicting results for the noise cu c dependence of G, 
both following Mott's work [TJ]. The question is whether to regard the noise as 
"low frequency driving" or as "temperature". On the one hand, on the basis of 
the Kubo formula, one expects a crossover from G ~ exp(— Ljl^ [where £oo is the 
localization length], to the noise dependent result G ~ lo 2 c \ log(w c )| d+1 , where d—1 for 
quasi one-dimensional (ID) ring. On the other hand, on the basis of the variable range 
hoping (VRH) picture, one expects G ~ exp(— (u)o/w c ) X ' d+1 ), where luq is a constant. 
Eventually [5] it has been realized that both the ballistic, the diffusive and the strong 
localization regimes should be handled on equal footing using SLRT. The Kubo theory 
applies in the LRT limit e — > 0, while in mesoscopic circumstances SLRT leads to a 
resistor network [13] "hopping" picture in energy space that generalizes the real space 
hopping picture of Refs.[l3 | ITB ] . 

2. Outline 

In this Communication we analyze, within the framework of SLRT, the dependence 
of the mesoscopic conductance of a quasi- ID ring on the strength of the disorder. We 
explain that for both weak and strong disorder the non-ergodicity of the quantum 
eigenstates implies having texture and sparsity in the perturbation matrix. Such 
features imply that the rate of energy absorption is suppressed enormously because 
the system cannot execute connected sequences of transitions. The implied deviations 
from the Kubo-Drude result are demonstrated numerically in Fig. 1. 

We introduce a novel random matrix theory (RMT) model, with either log- 
box or log-normal distributed elements, that captures the essential features of the 
perturbation matrix. A generalized resistor network analysis for the EMF-driven 
transitions in energy space leads to a generalized VRH picture in the strong disorder 
limit, while handling on equal footing the opposing limit of very weak disorder. 

In the first part of this communication (sections 1-5) we provide the essential 
details on the model and on the LRT/SLRT calculation, leading to the numerical 
results in section 6. The second part of this communication leads to the RMT modeling 
and to the implied generalized VRH picture. 

| The term "semi-linear response" to describe the outcome of the theory of Ref . [9] has been coined 
in a subsequent work where it has been applied to the analysis of the absorption of low frequency 
radiation by metallic grains. 
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3. Modeling 

We consider disorder quasi- ID ring geometry. The amount of disorder is traditionally 
characterized by the mean free path I. The semiclassical theory of the conductance 
leads to the Drude formula^ 

e 2 I 

G Drudo = — M - (2) 
znn L 

where L is the length of the ring, and M. is the number of open modes (proportional 
to its cross section). For the numerical calculations we have used the Anderson tight 
binding model, where the lattice is of size L x M with M <C L. The longitudinal 
hopping amplitude per unit time is cn = 1, while in the transverse direction it is 
numerically convenient to have c± < 1 , so as to have in the middle of the band a finite 
energy window with M. = M open modes. 

The random on-site potential in the Anderson tight binding model is given by a 
box distribution of width W. The density of states at the Fermi energy, and the mean 
free path in the Born approximation, are written as 

irnv F A V W J 

The implied definition of v F leads to its identification as the Fermi velocity in 
the absence of disorder (disregarding a prefactor of order unity). The Anderson 
localization length for L = oo sample would be £oo = M.I. Accordingly, for finite 
sample, depending on the strength of the disorder, one distinguishes between the 
ballistic regime (L <C i) the diffusive regime « L « i^) and the Anderson strong 
localization regime (L ^> 



4. The LRT calculation 



The Fluctuation-Dissipation version of the Kubo formula expresses the conductance 
as an integral over the velocity-velocity correlation function: 

/e\ 2 1 f°° 

G LRT = g F (- J X - J (v(t)v(0))dt (4) 

where v is the velocity in the longitudinal direction. The Drude formula Eq.@ is 
based on the simplest classical approximation: 

(«(<M0)> « « 2 exp[-2(^)|i|] (5) 

Our objective is to find the conductance of the closed ring in circumstances such that 
the motion inside the ring is coherent (quantum interferences are not ignored). The 
calculation involves the quantum version of (v(t)v(0)) which can be obtained as the 
Fourier transform of the spectral function 

C{w) = ^^| Wnm | 2 2^ r (^-(£; m -^„)) (6) 

nm 

% Optionally if the ring is characterized by its transmission g then l/L is replaced by g/(l—g). 
See Ref. [9] I10| for details. As could be expected the result is in agreement with the Landauer 
theory 1161 provided j C 1. Indeed the Landauer formula can be obtained from Kubo formula, 
using a semiclassical evaluation of velocity-velocity correlation function, whenever the contribution 
of trajectories with non-zero winding number can be neglected. 
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where N is the size of the energy window of interest. This spectral function can 
be re- interpreted as describing the band profile of the perturbation matrix {v nm }. 
In particular the calculation of C(0) = 27Tfj F ((|7; rim | 2 )) LBT involves a simple algebraic 
average over the near diagonal matrix elements at the energy range of interest. Using 
this notation the formula for the Kubo conductance takes the form 

Glrt = 7T qI ((|IW| 2 )) LRT (7) 

The 0(A/oj c ) weak localization corrections to the Drudc formula Eq.([2]) are 
determined by the interplay of the broadened delta function in Eq.© with the level 
statistics. See Ref.(3]. Equivalently we may say that the algebraic average ((• • -)) LRT 
has some weak sensitivity to the off-diagonal range of the averaging. 



5. The SLRT calculation 



As in the standard derivation of the Kubo formula, also within the framework of 
SLRT, the leading mechanism for absorption is assumed to be Fermi Golden Rule 
(FGR) transitions. These are proportional to the squared matrix elements |f nm | 2 of 
the velocity operator. Still, the theory of [9] does not lead to the Kubo formula. This is 
because the rate of absorption depends crucially on the possibility to make connected 
sequences of transitions. It is implied that both the texture and the sparsity of the 
|w nm | 2 matrix play a major role in the calculation of G. SLRT leads to a formula for G 
that can be cast into the form of Eq.fJTJ), provided the definition of ((...)) is modified. 
Following [HI [10] we regard the energy levels as the nodes of a resistor network. We 
define 

gnm = 2g- 3 J Vnm j\ 2 S r (E m -E n ) (8) 

Then it is argued that ((|vnm| 2 }}sLRT is the inverse resistivity of the network. It 
is a simple exercise to verify that if all the matrix elements are the same, say 
Km| 2 = o- 2 , then {{\v nm \ 2 }} 

slrt — G too. But if the matrix has texture or sparsity 

then {{\v nm \ 2 )} 

SLRT < «Kmi 2 » LRT ■ 



6. Numerical results 



It is natural to define the scaled conductance of the ring as follows: 

6 - wi§m - 2M x 4 <<K ™ |2>> (9) 

This would be the average transmission per channel, if we had considered open 
(Landauer) geometry. But for a closed (ring) geometry G is determined by the 
appropriate "averaging" procedure ((...)) LRT t SLRT . In the SLRT case the "averaging" 
is in fact a resistor network calculation. If all the near diagonal elements are 
comparable in size then SLRT will give essentially the same result as LRT. More 
generally ((|tWi| 2 ))sLRT is typically bounded from above by the algebraic average, and 
bounded from below by the harmonic average. The latter is defined as ((X)) b = 
and reflects "addition of resistors in series". If the distribution is not 
too stretched then the median, or the geometric average, or the mixed average of 
Ref.[9] might provide a good approximation. But in general a proper resistor network 
calculation is required. The resistor network calculation is sensitive to the texture 
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and to the sparsity of the perturbation matrix. By texture we mean that the gray- 
level image of the v nm matrix appears to be scarred by structures, rather than being 
homogeneous. Looking on the images of the v nm matrices for various values of W, one 
realizes that both texture and sparsity emerge in the ballistic case, while in the strong 
localization case one observes only sparsity. We further expand on the quantitative 
characterization below. 

In Fig.l we plot the Drude conductance G Drudo of Eq.Q, and the Kubo 
conductance G L r T , together with the mesoscopic conductance G SL rt versus W. We 
see that outside of the diffusive regime, for both weak and strong disorder, the SLRT 
result is extremely small compared with the LRT expectation. This generic behavior is 
related to the sparsity and to the texture of the perturbation matrix, which is implied 
by the statistical properties of the eigenstates. The statistical analysis is carried out 
in Fig. 2, while the RMT perspective is tested in Fig. 3. The content of Figs. 2-3 is 
further discussed in the following sections. 

In order to determine numerically whether the texture is of any importance we 
simply permute randomly the elements of the v nrn matrix along the diagonals, and 
re-calculate G (see the "untextured" data points in Fig. 3). Obviously, by definition 
G LRT is not affected by this numerical maneuver. But it turns out that Gslrt is 
somewhat affected by this procedure in the ballistic regime, but still the qualitative 
results come out the same. Accordingly we deduce that the main issue is the sparsity, 
and concentrate below on the RMT modeling of this feature. 

In the remainder of this paper, we pave an analytical approach to the calculation 
of the conductance, which will allow us to shed some light on these numerical findings. 
First we discuss the familiar diffusive regime where both LRT and SLRT should 
be in agreement with the Drude approximation (the latter should become a good 
approximation for a big sample). Then we discuss the departure of SLRT from LRT 
outside of the diffusive regime, which reflects the sparsity of the v nm matrix due to 
the non-ergodicity of the eigenfunctions. 

7. The Random Wave conjecture 

In the diffusive regime Mott has argued that the eigenstates of the Hamiltonian matrix 
are ergodic in position space, and look like random waves. Using this assumption one 
can reconstruct the Drude result. Following Mot10we assume that t is the correlation 
scale of any typical eigenfunction ty{x,y). The basic assumption of Mott is that the 
eigenstates are locally similar to free waves. The total volume L d is divided into 
domains of size £ d . Hence we have (L/l) d such domains. Given a domain, the condition 
to have non- vanishing overlap upon integration is \q n — q m \t < 2ir, where q is the local 
wavenumber within this domaim. The probability that q n would coincide with q m is 
1/ (kE() d ~ 1 ■ The contributions of the non-zero overlaps add with random signs hence 



where assuming ergodicity fy 2 ps 1/L d . From here we get G ~ IjL leading to the 
Drude result. We discuss the limited validity of this result in the next section. 

+ The original argument by Mott is somewhat vague. We thank Holger Schanz for helpful 
communication concerning a crucial step in the derivation. 




(10) 
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8. The non ergodicity issue 

It is clear that Mott's derivation of the Drude formula on the basis of LRT and 
the Random Wave conjecture becomes non-applicable if the eigenfunctions are non- 
ergodic. This is indeed the case for both weak and strong disorder: a typical 
eigenfunction does not fill the whole accessible phase space. In the ballistic regime a 
typical eigenfunction is not ergodic over the open modes in momentum space, while 
in the strong localization regime it is not ergodic over the ring in real space [17] . Fig. 2 
demonstrate this point by plotting the participation ratio as a function of the disorder 
strength. Lack of quantum ergodicity for either weak or strong disorder implies that 
the perturbation matrix v nm is very textured and/or sparse. For the following analysis 
a precise mathematical definition of sparsity is required. In the next sections we shall 
provide such a definition, but for this purpose we have to shed some light on the size 
distribution of the matrix elements. 

In the strong disorder regime the observed "sparsity" is very simple for 
understanding: Eigenstates that are close in energy are typically distant in real space, 
and therefore have very small overlap. The "big" matrix elements are contributed 
by eigenstates that dwell in the same region in real space, and hence sparse in 
energy space. What we are going to call in the next sections "sparsity", is merely 
a reflection of the associated log-box size distribution of the matrix elements (see 
Fig. 2b). The log-box distribution is deduced by a straightforward extension of the 
above argument. A generic eignefunction in the localized regime has an exponential 
shape il>(r) ~ exp(— \r — ro\/t oa ) which is characterized by the localization length 
Consequently a typical matrix element of {|w n m| 2 } has the magnitude 



where x £ [0, L/2] has a uniform distribution. The prefactor is most easily derived 
from the requirement of having (X) w (£/A4L)v% in agreement with the semiclassical 
result. The latter is deduced from the Fourier transform of the velocity- velocity 
correlation function Eq.([S]). 

In the weak disorder regime the explanation of the observed "sparsity" and 
textures requires some more effort. For the purpose of this Communication we shall 
be satisfied with a qualitative explanation: If the disorder W were zero, then the mode 
index (call it n v ) would become a good quantum number. This means that states that 
are close in energy are not coupled (because they have different n y ). Once W becomes 
non-zero (but still small) the mixing is described by Wigner Lorentzians (much the 
same as in the toy model of Ref. [9]). Then the ratio between small and large couplings 
is determined by the different degree of mixing of close versus far modes. Consequently 
one observes a wide (but not stretched) distribution for the \og(X) values (see Fig. 2b). 

9. RMT modeling, beyond the Gaussian assumption 

It was the idea of Wigner [T5] to model the perturbation matrix of a complex system as 
a (banded) random matrix. Later it has been conjectured by Bohigas [19) that similar 
modeling may apply to chaotic systems. For many purposes it is convenient to assume 
infinite bandwidth with elements that are taken from a Gaussian distribution, leading 
to the standard Gaussian Orthogonal or Unitary ensembles (GOE/GUE). But there 
are obviously physical circumstance in which it is essential to go beyond the Gaussian 



X 
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assumption, and to take into account the implications of having finite bandwidth 
and/or non-Gaussian distribution of elements [30] and/or sparsity [H] and/or texture. 

It should be clear that the default assumption of having a Gaussian distribution 
of in-band matrix elements is legitimate on practical grounds as long as the matrix 
elements have comparable size in absolute value. But if the eigenfunctions are non- 
ergodic this assumption becomes problematic, because the elements (in absolute value) 
might have a wide distribution over many decades in log scale. In such case different 
type of averages may differ by orders of magnitude. 

In the following we regard {|unm| 2 } as a random matrix of non-negative 
numbers {X}. In general it might be a banded matrix. If the standard Gaussian 
assumption applies then the in-band elements of {X} are characterized by the Proter- 
Thomas distribution. But we are interested in physical circumstances in which many 
of the in-band elements are vanishingly small. We define this feature as "sparsity" . 
In the next section we define p as the fraction of elements that are larger than the 
average. If we have p <C 1 then we say that the matrix is "sparse" . We further discuss 
the definition of p in the next section. 

10. Characterization of sparsity 

For an artificially generated sparse random matrix {X} of non-negative elements, 
one defines p as the fraction of non-zero elements. Such a definition assumes a 
bimodal distribution. But in general realistic circumstances we do not have a bimodal 
distribution. Rather for strong disorder we already had explained that the distribution 
of the matrix elements {|v nm | 2 } is log-box. Contemplating a bit on this issue 
one concludes that the physically generalized definition of the sparsity measure is 
p = F((X)) where F(X) is the probability to find a value larger than X. We regard 
a matrix as sparse if p <C 1. Given that ln(X) is uniformly distributed between 
ln(Xo) and ln(Xi) we define p= (ln(Xi/ Xq))^ 1 , and find assuming Xq<^Xi that 
p w —p\np, and (X) fa pX\. Hence for log-box distribution (X) ~ pX\, as expected 
from the standard bimodal case. 

In Fig. 3 we re-do the calculation of the conductance with artificial matrices with 
the same sparsity, i.e. log-box distributed elements with the same p. We observe 
qualitative agreement for strong disorder. In the other extreme limit of weak disorder 
there is no agreement, because we have to use a different distribution for the matrix 
elements: It turns out that also in the ballistic regime log(X) has a wide distribution, 
but it is not stretched as in the case of a log-box distribution (see Section 8). In practice 
we can describe the X distribution of the matrix elements in the ballistic limit as log- 
normafl Once we use the appropriate distribution we get a reasonable qualitative 
agreement. We emphasize that in both cases, of either weak or strong disorder, 
there is besides the algebraic average (X) only one additional fitting parameter that 
characterizes the distribution and hence determines the "sparsity" . We could of course 
have generated RMT matrices using the actual distribution (Fig. 2), but then we would 
merely re-generate the untextured data points. 

* The default fitting of the log(X) distribution to a Gaussian line shape is merely a practical issue. 
The "RMT ideology" is to see whether a "minimum information" ensemble of random matrices can 
be used in order to derive reasonable estimates. If we want to further improve our estimates in 
the ballistic regime it is essential to take into account the texture and not only the deviation from 
log-normal distribution (see Section 6). 
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Given a hopping range \E m — E n \ < lu we can look for the typical matrix element 
X for connected sequences of such transitions, which we find by solving the equation 

(s) F « ~ 1 W) 

In particular for strong disorder we get 

x « ^^/M (13) 



G 



where = (L/l^A is the local level spacing between eigenstates that are localized 
in the same region. The same procedure can be applied also in the ballistic regime 
leading to a simpler variation of l|13p where the dependence of X on uj predominantly 
reflects the band profile: It follows from the discussion after Eq.® that v nm is a 
banded matrix, with a Lorentzian bandprofile whose width ~ v F /£ becomes narrower 
as the disorder is decreased. 

11. Generalized Kubo formula 

The definition of the bandprofile reflects the variation of (X) with to. In complete 
analogy we define an effective bandprofile that reflects the variation of X with to. 
Namely: 

C qm _ LRT H = 2irg F (X) (14) 

CWslrtM = 2irg F X (15) 

The spectral function of Eq.© is a smeared version of the "bare" spectral function: 
it is obtained by a convolution C qm _ LRT (cj) *5r(uj). Consequently we get 

5 (J) 2 Q F j ' CMSrWfo (16) 

where the appropriate LRT/SLRT spectral function should be used. This way of 
writing allows to obtain an approximation for the mesoscopic conductance using a 
Kubo-like calculation: it is just re-writing of the Kubo formula in the LRT case, while 
being a generalized VRH approximation in the SLRT case. 

For strong disorder the above generalized VRH approximation gives an integral 
over exp(— |w|/w c ) exp(— Ae/\u>\), which is a product of two competing factors: the first 
has to do with the noise/temperature and the second has to do with the couplings. 
This integral is handled using the usual VRH phenomenology: the result is determined 
by the maximum of its integrand, which requires to optimize the range u> of the 
transition. In the weak disorder regime the VRH integral is not the same as in the 
strong disorder case, because a log-normal rather than log-box distribution is involved. 
We have verified that the generalized VRH integral gives a qualitatively reasonable 
approximation to the actual resistor network calculation in both cases. In any case 
one should keep in mind the well known reservations that apply to such "mean field" 
approach [T4l [T5] . 

12. Summary 

Within SLRT it is assumed that the transitions between levels are given by the Fermi- 
golden- rule (FGR), but a resistor network analogy is used in order to calculate the 
energy absorption rate. The calculation generalizes the variable range hopping picture 
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and treats on equal footing the weak and strong disorder regimes. The essential physics 
is captured by RMT, provided the perturbation matrix is regarded as a member of 
the appropriate Gaussian / Log-Normal / Log-Box ensemble. 

The prevailing results in the literature regarding the conductance of small closed 
metallic rings (for a review see Ref. [4]) concern mainly the diffusive regime, where 
in leading order the conductance is given by the Drude formula, and SLRT does not 
differ much from LRT. In the present communication multi-mode rings in the non- 
diffusive regime are considered seriously for the first time. Then it become essential 
to define the precise assumptions regarding the environment and the driving. It is 
important to realize that both LRT and SLRT assume Markovian FGR transitions. 
This is a very realistic assumption that can be justified rigorously if one assumes noisy 
driving (as in our exposition) or else it is implied by having a noisy environmenl|jJ. 
Accordingly it should be clear that LRT and SLRT both share the same small 
parameter as in the FGR picture, which is the ratio between the rate of the driven 
transitions and the smallest relevant energy scale that characterizes the band-profile 
(bandwidth / sparsity/texture) . 

There is only one assumption that distinguishes the SLRT (mesoscopic) 
circumstances from LRT (Kubo) circumstances. This is related to the implicit role of 
the environmentally induced relaxation process in the determination the steady state 
of the system. Within SLRT one assumes that the FGR rate of the driven transitions 
(iu FGR oc e 2 g nm ) is larger compared with the relaxation rate (7 r]x ). The inelastic 
relaxation effect can be incorporated into the SLRT framework by considering a non- 
symmetric g rml as implied (say) by detailed balance considerations. If the relaxation 
is the predominant effect (w FGR < 7) then we are back in the LRT regime (23] where 
the Kubo-Drude result applies [S]. 

One can wonder what happens if the FGR assumption of LRT/SLRT breaks down. 
Not much is known [24] . Ref. [25] has attempted to go beyond the FGR approximation 
using the Keldish formalism, and has recovered a Marokovian picture that leads to 
a Kubo-like result for the conductance. If the Keldish Markovian picture could be 
established beyond the diffusive regime , it would be possible to extend SLRT into 
the non-linear regime. 
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Fig.l: Plot of the scaled conductance G versus W using either the LRT (Kubo) or the SLRT 
(mesoscopic) recipe, and compared with Drude formula based estimate. The calculation has been 
carried out for a tight binding Anderson model of size 500 X 10, transverse hopping amplitude c=0.9 
and low driving frequency oj c /A = 7. The SLRT result departs from the LRT result for both weak 
disorder (ballistic regime) and strong disorder (strong localization regime). 
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Fig. 2: The ergodicity of the eigenstates is characterized by the participation ratio PR = p 2 ] -1 , 
which is calculated (left panel) in various representations: in position space p rx ,r y = {r x > °r y \ 2 , 
in position-mode space p T:c ,ky = \( r x, k y \^)\ 2 , and in mode space p^ y = \{r x , ky\¥>)\ 3 , where 

k y = [-k/(M+1)\ X integer. The cumulative distribution F(X) of the in-band matrix elements 
(right panel) exhibits a log-box distribution in the strong localization regime. Points in the interval 
X > ((X))slrt , corresponding to non-negligible values, are connected by a thicker line. The extracted 
sparsity measure (left panel) is p = F((X)). 
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Fig. 3: The ratio Gslrt/Glrt versus the inverse of GDrudc = £/L based on the numerics of Fig.l, 
and compared with artificial RMT modeling using "sparse" matrices formed of log-normal or log-box 
distributed elements. We also compare the actual results with "untextured" results as explained in 
the text. For weak disorder the agreement is only qualitative indicating that the texture becomes 
important. 



